Search for the decay n^uu in the momentum region < 195 MeV/c 
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We have searched for the decay K'^ —* -n^vv in the kinematic region with pion momentum 
below the — > tt^tt" peak. One event was observed, consistent with the background estimate of 
0.73 ± 0.18. This implies an upper limit on B{K+ -k+vv) < 4.2 x 10"^ (90% C.L.), consistent 
with the recently measured branching ratio of (l-57/Qg2) x 10"^^°, obtained using the standard 
model spectrum and the kinematic region above the K'^ tt'^tt'^ peak. The same data were used 
to search for —* n'^ X'^, where is a weakly interacting neutral particle or system of particles 
with 150 < Mxo < 250 MeV/c^ 

PACS numbers: 13.20.Eb, 12.15.Hh, 14.80.Mz 



In a recent paper we reported the branching ratio for 
the rare decay K+ tt+w to be (1.57tJ;^^) x 10"^" 
based on the observation of two events in the phase space 
region > 211 MeV/c |l[. This decay is sensitive to 
the coupling of top to down quarks, Vtd, in the Cabibbo- 
Kobayashi-Maskawa mixing matrix. The standard model 
(SM) predicted branching ratio, B{K^ tt^vv), is 
(0.75 ± 0.29) X lO^^"!^. Loop diagrams involving new 
heavy particles in extensions of the SM can interfere with 
SM diagrams and alter the decay rate, and also the kine- 
matic spectrumQ. Exotic scenarios such — > tt+X^ 
where X'^ is a hypothetical stable weakly interacting par- 
ticle or system of particles have also been suggested |^, |j . 
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It is therefore important to obtain higher statistics for 
this decay and to extend the measurement to other re- 
gions of phase space. The results in ||] are from analysis 
of data with the tt"*" momentum above the —> tt+tt" 
{KTr2) peak (Region 1). The 7r+ from Kt^2 decay has a 
kinetic energy (E), momentum (P), and range (R) in 
plastic scintillator of 108 MeV, 205 MeV/c, and 30 cm, 
respectively. In this letter we report the analysis of data 
below the ifT2 peak (Region 2) obtained from Experi- 
ment E787§, 0, |, |, |l^, y, |l| at the Alternating Gra- 
dient Synchrotron (AGS) of Brookhaven National Lab- 
oratory. The previous limit using Region 2, 1.7 x 10~^ 
(90% C.L.), was obtained from an earlier version of the 
E787 detector (l|l. 

The signature for tt'^vv in the E787 experi- 

ment is a single K'^ stopping in a target (TG), decay- 
ing to a single ■k'^ with no other accompanying pho- 
tons or charged particles. In Region 1, the major back- 
grounds were found to be the two body decays Kt;2 
and —>■ ijl^v^ (^^^2), scattered beam pions, and 
charge exchange (CEX) reactions resulting in decays 
TT~^l~Di, where Z = e or /i. Region 2 has larger po- 
tential acceptance than Region 1 because the phase space 
is more than twice as large and the loss of pions due to 
nuclear interactions in the detector is smaller at the lower 
pion energies. However, there are additional sources of 
background for Region 2. These include in which 
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the 7r+ loses energy by scattering in the material of the 
detector (primarily in the TG), tt^tt'^j (-^^7727)1 

K+ ^ fi+iyj {K^2^), K+ ^+i.^o (i^^a), and K+ ^ 
(ife4) decays in which both the tt^ and the 



e"*" are invisible because of absorption. 

The data were obtained with a flux of 6 x 10® kaons 
per 1.6 sec spill at 730 MeV/c (with 24% pion contam- 
ination) entering the apparatus. The kaons were identi- 
fied by a Cerenkov detector; two multi- wire-proportional- 
chambers were used to determine that there was only one 
entering particle. After slowing in a BeO degrader the 
kaons traversed a 10-cm-thick lead-glass detector read 
out by 16 fine-mesh photomultiplier tubes (PMT) and 
a scintillating target hodoscope (TH) placed before the 
TG. The lead-glass detector was designed to be insensi- 
tive to kaons and detect electromagnetic showers origi- 
nating from kaon decays in the TG. The TH was used to 
verify that there was only one kaon as well as determine 
the position, time, and energy loss of the kaon before 
it entered and stopped in the TG. The TG consisted of 
413 5.0-mm-square, 3.1-m-long plastic scintillating fibers, 
each connected to a PMT. The fibers were packed axi- 
ally to form a cylinder of '--^12 cm diameter. Gaps in the 
outer edges of the TG were filled with smaller fibers which 
were connected to PMTs in groups. The PMTs were read 
out by ADCs, TDCs, and 500 MHz transient digitizers 
based on GaAs charge-coupled devices (CCDs)[0. Pho- 
tons were detected in a hermetic calorimeter mainly con- 
sisting of a 14-radiation length thick barrel detector made 
of lead/scintillator sandwich and 13.5-radiation length 
thick endcaps of undoped Csl crystals ||] . The rest of the 
detector consisted of a central drift chamber (UTC)j9|, 
and a cylindrical range stack (RS) of 21 layers of plastic 
scintillator with two layers of embedded tracking cham- 
bers, all within a 1-T solenoidal magnetic field. The TG, 
UTC, and RS allowed the measurement of the P, i?, and 
E of the charged decay products. The tt ^ ^ e decay 
sequence from pions that came to rest in the RS was ob- 
served using another set of 500 MHz transient digitizers 
(TD)0. 

The data reduction and offline analysis for Region 2 
was similar to the analysis of Region 1 |l^, [ij, |lj] , al- 
though the final cuts to enhance signal and suppress 
background to less than one event were different. Here 
we will emphasize the key instrumentation and analy- 
sis tools used to suppress the background in Region 2. 
The TG, CCDs, and the photon veto system were the 
important elements for Region 2 analysis. A multilevel 
trigger selected events by requiring an identified K'^ to 
stop in the TG, followed, after a delay of at least 1.5 ns, 
by a single charged particle track that traversed TG and 
RS with a hit-pattern consistent with the expectation for 
TT'^vD. Events with photons were suppressed by 
vetos on the barrel and endcap detectors. In the of- 
fline analysis, the single charged particle was required to 
be identified as a tt"*" with P, R, and E consistent for 
a TT"*", and the TD pulse information consistent with the 
decay sequence tt — > /i — > e in the last RS counter on 
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FIG. 1: Range (cm in plastic scintillator) and kinetic energy 
(MeV) of events remaining after all cuts except the momen- 
tum cut (top), and including the momentum cut (bottom). 
The dark points represent the data. The simulated distribu- 
tion of expected events from K'^ -k^vu is indicated by the 
light dots. The group of events around 108 MeV is due to the 
K-n2 background. The events at higher energy are due to _K'p2 
and Kfj_2-f background. All events except for the one in the 
signal box are eliminated by the 140 < P < 195 MeV/c cut 
on momentum. 

the pion trajectory. The signal region was defined by the 
intervals 140 < P < 195 MeV/c, 12 < i? < 27 cm, and 
60 < £; < 95 MeV. 

The background was found to be dominated by Kj^2 
events in which the pion had a nuclear interaction near 
the kaon decay vertex, most probably on a carbon nu- 
cleus in the TG plastic scintillator. This scatter left the 
pion with reduced kinetic energy, putting it in Region 
2. We suppressed this background by removing events 
in which the pion track had a scattering signature in the 
TG. These signatures included kinks, tracks that did not 
point back to the vertex fiber in which the kaon decayed, 
or energy deposits inconsistent with the ionization en- 
ergy loss for a pion of the measured momentum. The 
remaining iir^2 background consisted of events in which 
the pion scattered in one of the fibers traversed by the 
kaon. The extra energy deposits from the pion scatters 
were obscured by the earlier large energy deposits of the 
kaon. For these events, we examined the pulse shapes 
recorded in the CCDs in each kaon fiber using a fit and 
eliminated events in which an overlapping second pulse, 
in time with the pion, was found to have energy larger 
than 1 MeV. To obtain sufficient separation of the 
and TT^ induced pulses in the CCDs we required a min- 
imum delay of 6 ns between the kaon and the pion. Pi- 
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K+ tt+tt" 


d 


0.630 ±0.170 


K+ -> 7r+7r% 


dm 


0.027 ± 0.004 




d 


0.007 ± 0.007 


Beam 


d 


0.033 ± 0.033 




dm 


0.026 ± 0.032 


CEX 


dm 


0.011 ±0.011 


Total 




0.734 ±0.177 



TABLE I: Estimated number of background events for Re- 
gion 2 of K'^ — > -K^vv data. The second column indicates the 
method of background determination: data alone (d), data 
combined with simulation (dm). The errors include statistics 
of the data and Monte Carlo as weU as systematic uncertain- 
ties. 

nally, additional Kt^2 rejection was obtained by removing 
events with photon interactions in detectors surrounding 
the kaon beam-hne; these cuts caused substantial (~ 42 
%) loss of efficiency because of accidental hits due to the 
high flux of particles. 

We formed multiple independent constraints on each 
source of background. These constraints were grouped 
in two independent sets of cuts, designed to have little 
correlation. One set of cuts was relaxed (or inverted) to 
enhance the background so that the other set could be 
evaluated to determine its power of rejection, as summa- 
rized below. The sum of the background due to K^2-y 
and K^'i was obtained by separately measuring the re- 
jection factors of the TD particle identification and kine- 
matic {R and P) particle identification. The background 
from beam pions was evaluated by separately measuring 
the rejections of Cerenkov, TH beam particle identifica- 
tion, and the delay time between pion and kaon. The 
dominant background from Kt^2 decay was measured by 
evaluating the rejection of the photon veto system on 
events tagged by scattering signatures in the TG and 
target CCDs. Similarly, the rejection of the target CCD 
cut was determined by using events that failed the photon 
veto criteria, ft should be noted that the Region 1 analy- 
sis measured photon veto rejection using the unscattered 
events in the momentum peak (205 MeV/c)[|l^. This 
method could not be used for Region 2 because the scat- 
tering in the TG spoiled the back-to-back correlation be- 
tween the detected 7r+ track in the RS and the undetected 
tt", leading to different photon veto rejection factors for 
scattered and unscattered Kt^2 background events. 

We employed Monte Carlo (MC) simulation to evalu- 
ate the backgrounds from Kf,^ and i4r^27 because these 
could not be distinguished, on the basis of the 7r+ track 
alone, from the much larger Kt^2 background with a 7r+ 
scatter. In the case of K^i, we first identified such events 
by looking for additional short tracks, indicative of tt^ 
or , coming from the kaon decay vertex in TG. We 
then used MC events to estimate the power of rejecting 
background events with charged particle hits unrelated 
to the K+ or the 7r+ in the TG. The MC simulation 
was performed using the previously measured decay dis- 
tribution of Kf-i decays For the final estimate (Ta- 
ble ^ we relied on the observed number of Ke4^ candi- 



Acceptance factors 




stop efficiency 


0.670 


Ji"*" decay after 6 ns 


0.591 


K'^ — > ■k'^vv phase space 


0.345 


Geometry 


0.317 


TT^ nucl. int. and decay in flight 


0.708 


Reconstruction efficiency 


0.957 


Other kinematic cuts 


0.686 


n — fj, — e decay chain 


0.545 


Beam and target analysis 


0.479 


CCD acceptance 


0.401 


Accidental loss 


0.363 


Total acceptance 


7.65 X 10"* 



TABLE IL Acceptance factors used in the measurement of 
-K^vv in Region 2. The "A'^ stop efficiency" is the 
fraction of kaons entering the TG that stopped. "Other kine- 
matic constraints" include particle identification cuts. 

date events, which were found to be consistent with the 
known branching ratio. For simulation of tt" absorption 
we used a measurement of ionization spectra in scintil- 
lator performed previously by our collaboration |lq| and 
corroborating information from other sources p^ . For 
the if,r27 background, a MC calculation provided the ra- 
tio between the observed number of events in the Kt^2 
peak and the expected background. The calculation was 
performed using the previously measured [^| decay rate 
and a calculation of the extra rejection, based on calibra- 
tion data, due to the radiative photon in isr^27 decays. 
Measurements of K'^ charge exchange reaction in the TG 
were performed. These measurements were used as input 
to the MC calculation of the CEX background. The final 
background estimates and associated errors in Table | 
include corrections for small correlations in the separate 
groups of cuts and cross contamination of background 
samples. 

The integrity of the background estimates was assured 
because the background cuts were defined using only one- 
third of the data, sampled uniformly from the entire set, 
without examining the events in the pre-determined sig- 
nal region. The cuts were then applied with no further 
changes to the remaining two-thirds of the data to ob- 
tain the numbers reported in Table |. The systematic 
error on the largest background, ^^^2, was estimated by 
measuring the rejection of photon veto cuts on many dif- 
ferent event ensembles, tagged in different ways for a TG 
scattering signature. The event ensembles were designed 
to have little contamination by other background sources 
such as if7r27 and iCe4- As a final check, each cut was re- 
laxed to admit background events in a predictable way. 
Examination of these background events, which are close 
to the signal region, provided no indication of background 
sources other than those in Table |. For example, the 
kaon decay time region between 2 to 6 ns, with accep- 
tance of 0.254 ± 0.004 (less than the naive expectation 
due to other lifetime-dependent cuts) relative to the sig- 
nal region, was examined. This region has a total es- 
timated background of 0.45 ± 0.14, dominated by Kt^2 
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FIG. 2: The 90% C.L. upper limit for B{K+ ~+ n+X^) as a 
function of Mxo, the mass of the recoihn g sy stem. The sohd 
(dashed) line is from this analysis (from 113]). The limit for 
Mxo < 140 MeV/c^ is derived from the result for Region 1. 
The observation of one event in Region 2 reported in this Let- 
ter causes a bump in the limit at 194 MeV/c^. Similarly the 2 
events, consistent with the observation of K'^ -k^uv above 
background, in Region 1 reported in |l| increase the limit at 
105 and 86 MeV/c^. We have also included the single event 
sensitivity as a function of Mxo (the dotted line), defined in 
the text, obtained by E787. 

decays due to the reduced background rejection from the 
CCDs. One event, consistent with the background esti- 
mate, was observed in this background region. 

After the background study, the signal region was 
examined, yielding one candidate event with P = 
180.7 MeV/c, R = 22.1 cm, and E = 86.3 MeV with 
a kaon decay time of 17.7 ns, consistent with the back- 
ground estimate of 0.73 ± 0.18. Fig. |l| shows the kine- 
matics of the remaining events before and after the cut 
on measured momentum, P. 

Using the total number of incident on TG for these 
data, 1.12 x 10^^, the acceptance reported in Table||, and 
the observation of one event in Region 2 we calculate the 
upper limit of B{K+ ir+vv) < 4.2 x 10"^ (90% C.L.) 



p0{ . This is consistent with the branching ratio reported 
from Region 1 and the SM decay spectrum combin- 
ing the measurements from Region 1 and Region 2 does 
not alter the branching ratio measurement significantly 
because it is dominated by the sensitivity of Region 1. 
However, for non-standard scalar and tensor interactions. 
Region 2 has larger acceptance than Region 1. We have 
combined the sensitivity of both regions to obtain the 
90% C.L. upper limits, 4.7 x 10"^ and 2.5 x 10^^ for 
scalar and tensor interactions, respectively. 



This measurement is also sensitive to 



where X° is a hypothetical stable weakly interacting par- 
ticle, or system of particles. Fig. || shows 90% C.L. up- 
per Hmits on B{K~^ — > 7r+X°) together with the pre- 
vious limit from |l^. The dotted line in Fig. || is the 
single event sensitivity defined as the inverse of the ac- 
ceptance for K'^ TT+X" multiplied by the total number 
of stopped kaons as a function of Mxo. 

In conclusion, the use of GaAs charged-coupled devices 
to record pulse shapes as well as highly efficient photon 
detection has allowed us to suppress background in Re- 
gion 2. This has resulted in new limits on the spectrum 
of the pion in the decay — > n'^h'iy as well as im- 
provement in the sensitivity to A'+ —>■ tt+A^ by a factor 
between 4 and 40 over the accessible mass range. The 
detailed enumeration of backgrounds in Region 2 will be 
important for new experiments that intend to precisely 
measure t^^vD with large statistics [^l], p2| . 
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